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ABSTRACT 

The  heat  transfer  characteristics  of  a  container  stored  rocket 
motor  have  been  investigated  using  analytical  and  experimental 
techniques.  Comparison  between  analytically  predicted  and  experi- 
mentally determined  values  of  temperature  are  within  the  estimated 
experimental  uncertainty  of  ±3°F.  The  results  of  the  analytical 
solution  may  be  used  to  predict  maximum  and  minimum  temperatures, 
thermal  time  lags,  and  temperature  gradients  throughout  the  rocket 
motor.  For  convenience,  maximum  temperature  and  time  lag  Information 
is  presented  in  nomograph  form  as  a  function  of  three  parameters: 
normalized  radial  location,  normalized  frequency,  and  normalized 
gap  conductance.  These  three  parameters  control  the  temperature 
behavior  of  the  motor.  It  is  proposed  that  the  nomographs  will  be 
a  useful  tool  for  thermally  optimizing  future  container  designs. 


This  task  was  supported  by:      Naval  Weapons  Center,  China  Lake, 

California,  Work  Request  No.  3-3006 
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I.  INTRODUCTION 

An  analytical  model  is  presented  for  predicting  the  daily  tempera- 
ture history  of  a  container  stored  rocket  motor.  The  model,  which  is 

* 
developed  using  the  method  of  complex  temperatures  [1,2]  ,  assumes 

that  heat  is  transferred  only  in  the  radial  direction  and  that  the 
container  surface  temperature  is  known  a  priori  and  varies  si nusoi dally 
with  time.  Comparison  of  theory  with  experiment  is  within  experimental 
uncertainty  (±  3°F)  when  temperature  is  interpreted  as  bulk  temperature. 
The  analytical  model  is  especially  useful  for  studying  geometrical  and 
thermophysical  property  effects  on  rocket  motor  temperature.  Such  para- 
meter studies  have  been  carried  out  and  the  results  are  presented  in 
nomograph  form. 

II.  BACKGROUND 

In  1959  the  Naval  Weapons  Center  located  at  China  Lake,  California, 
recognized  the  need  for  a  concerted  attack  on  the  problem  of  thermal 
criteria  assignment  for  new  weapon  systems.  In  1963  a  task  force  was 
established  to  study  the  complete  environmental  criteria  determination 
problem.  The  key  to  this  problem  seemed  to  be  the  thermal  area  in  the 
storage  and  transportation  events  of  any  item.  It  was  realized  that 
transportation  was  a  short  term  situation  compared  to  the  storage 
situation.  Therefore,  the  major  portion  of  the  life  of  an  item  must 
be  in  storage.  There  are  three  types  of  storage;  covered,  igloo,  and 
dump.  The  dump  storage  situation  leads  to  the  more  extreme  thermal 
exposure  situations. 


* 
Numbers  in  brackets  designate  references  at  end  of  report. 


As  data  were  not  available  for  the  dump  storage  situation,  instru- 
mented storage  dumps  were  created  at  representative  places  on  a 
worldwide  basis  so  that  statistical  data  could  be  obtained  on  a  variety 
of  ordnance.  The  first  site  was  at  China  Lake,  California,  in  the 
middle  of  the  Mojave  Desert.  This  site  now  has  the  capability  of 
returning  about  250  channels  of  information  on  a  continuous  temperature- 
time  basis  (Figure  1).  Other  arctic  and  tropical  sites  were  set  up  to 
study  extreme  conditions. 

In  the  experimental  dump  storage  situation  the  ordnance  is  exposed 
singly,  directly  situated  on  the  ground,  with  the  long  axis  aligned  in 
the  north-south  direction  to  allow  maximum  normal  exposure  of  the 
container  surface  to  the  sun's  rays.  In  actual  practice,  ordnance  is 
usually  stacked  and  oriented  in  other  than  a  north-south  direction, 
thereby  avoiding  the  extreme  temperature  condition.  Ordnance  sitting 
on  the  ground  receives  terrestial  radiation,  cannot  quickly  give  off 
heat  by  conduction  to  the  soil,  and  is  not  as  apt  to  be  cooled  by  the 
prevailing  breeze;  therefore,  extreme  temperatures  result. 

The  most  important  source  of  heat  to  the  ordnance  is  the  direct 
radiation  from  the  sun,  with  terrestial  radiation  being  of  secondary 
importance.  For  extreme  conditions  to  occur,  the  wind  must  be  calm 
(less  than  5  knots),  the  sky  clear,  and  the  outside  air  temperature 
high.  After  sunrise  the  ordnance  skin  temperature  rises  much  more 
rapidly  than  the  ambient  air  temperature;  therefore,  the  surrounding 
air  cools  the  ordnance,  rather  than  heats  it. 

The  rocket  motors  used  for  the  tests  were  military  surplus.  Even 
though  the  material  had  served  its  intended  in-Fleet  purpose,  it  was 
still  representative  of  new  hardware,  when  viewed  in  a  thermodynamic 
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context.  When  inert  rocket  motors  were  available,  they  were  used  intact; 
however,  in  most  cases,  once-fired  hardware  was  used.  Thoroughly  dried 
desert  blown  sand,  being  similar  in  thermal  properties  to  most  propel- 
lants,  was  used  to  backfill  empty  rocket  motors.  It  was  assumed  that 
the  thermal  response  of  the  sand  filled  motors  was  essentially  the  same 
as  actual  propel  1  ant  filled  motors. 

III.  EXPERIMENTAL  PROCEDURE 

Although  the  Naval  Weapons  Center,  China  Lake,  had  accumulated  vast 
amounts  of  data  in  the  past,  it  was  decided  to  instrument  a  rocket  motor 
storage  container  system  especially  for  this  project.  This  would  allow 
base  data  to  be  taken  exactly  when  and  where  it  was  required.  It  also 
allowed  variations  in  the  system  without  interfering  with  one  of  China 
Lake's  ongoing  projects.  An  ASROC  (antisubmarine  rocket)  system  was 
chosen  for  this  study.  The  outer  storage  container  was  75  inches  long 
with  an  inner  diameter  of  18  inches  and  a  wall  thickness  of  1/16  inch. 
The  rocket  motor  was  57  inches  long  with  an  outside  diameter  of  12  inches 
and  a  wall  thickness  of  1/4  inch.  Both  the  container  and  motor  were 
made  of  steel. 

The  rocket  motor  storage  container  system  was  instrumented  with 
20  gage  copper-constantan  insulated  thermocouple  wire  with  an  ISA 
calibration  of  ±  1-1 /2°F  over  the  range  -75  to  +200°F.  Twenty-one 
thermocouples  were  originally  placed  on  the  system.  Positions  are 
indicated  in  Figures  2  and  3.  The  ambient  air  temperature  was  measured 
with  thermocouple  number  19  which  was  located  in  a  Stevenson  shelter 
about  60  feet  away  from  the  system  (Figure  4). 


Figure  2.   Thermocouple  Locations  on  Experimental  System. 


Five  thermocouples  were  located  under  the  section  painted 
with  the  liquid  crystals.   Their  locations  corresponding  to 
the  ones  shown  above  are:  #14=  #1,  #15=  #2,  #16=  #8,  #17=  #9, 
and  #18=  #3  (See  Figure  3) . 
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Figure  4.   Stevenson  Shelter. 


The  thermocouples  were  mounted  intrinsically  on  the  motor  and 
storage  container.  Two  small  holes  were  drilled  approximately  1/8 
inch  apart  in  the  metal  and  the  individual  wires  were  inserted  in  the 
holes.  The  metal  was  then  hammered  around  the  wires  until  a  snug  fit 
was  obtained.  Bead  thermocouples  were  mounted  at  the  center  of  the 
motor  and  in  the  air  gap.  The  thermocouples  located  at  the  center  of 
the  motor  were  supported  by  small  pieces  of  wood  several  inches  from 
the  bead.  The  use  of  these  supports  was  necessary  to  keep  the  thermo- 
couples in  position  when  the  motor  was  filled  with  sand.  After  all 
the  thermocouples  on  the  rocket  motor  were  in  place,  the  rocket  motor 
was  filled  with  dry  desert  blown  sand.  The  wires  from  the  two  thermo- 
couples located  in  the  center  of  the  motor  were  led  out  a  hole  in  the 
end  cap.  To  avoid  settling  of  the  sand  after  the  motor  was  in  place 
on  the  site,  with  a  resulting  air  gap  being  formed  between  the  sand 
and  the  motor  skin,  the  sand  was  compacted  by  striking  the  sides  of 
the  motor  with  small  sledge  hammers  and  then  adding  additional  sand 
through  the  hole  in  the  end  cap.  This  was  continued  until  the  sand 
was  tightly  packed.  The  hole  in  the  end  cap  was  then  sealed. 

The  rocket  motor  was  carefully  placed  in  its  storage  container 
(Figure  5)  which  had  previously  been  instrumented  with  thermocouples. 
The  thermocouples  in  the  air  gap  were  mounted  by  affixing  the  lead  wire 
to  the  rocket  motor  at  the  desired  position.  A  90  degree  bend  was  then 
put  in  the  wire  so  that  it  placed  the  bead  of  the  thermocouple  approxi- 
mately 1.5  inches  into  the  air  gap.  Neither  the  thermocouples  in  the 
center  of  the  motor  nor  those  in  the  air  gap  could  be  considered 
accurately  positioned.  Every  effort  was  made,  however,  to  minimize 
positioning  errors.  All  thermocouple  wires  were  located  inside  the 
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storage  container  and  were  led  through  a  hole  in  one  end.  This  hole 
was  then  sealed.  The  two  halves  of  the  storage  container  were  then 
bolted  shut. 

The  outer  surface  of  the  rocket  motor  and  the  inner  and  outer 
surfaces  of  the  storage  container  were  all  painted  various  shades  of 
haze  gray.  Weathering  had  caused  the  painted  surfaces  to  appear  fairly 
rough.  This  is  typical  of  the  conditions  of  a  storage  dump.  From  the 
condition  of  the  surfaces,  it  was  estimated  that  the  emissivity  was 
approximately  0.9. 

Prior  to  loading  the  rocket  motor  into  the  storage  container, 
liquid  crystals  were  applied  to  part  of  the  storage  container  surface 
in  order  to  obtain  a  thermal  mapping  of  the  surface  temperature  at  any 
instant  of  time.  Liquid  crystals  are  temperature  sensitive  materials 
that  produce  immediate  thermal  images  in  a  pattern  of  colors  which 
respond  rapidly  to  minute  changes  in  substrate  surface  temperatures. 
Prior  to  applying  the  crystals,  a  15  inch  strip  of  the  storage  container, 
20  inches  from  one  end,  was  sprayed  with  two  coats  of  Testors  Spray  PI  a 
Enamel  No.  1249,  Flat  Black  as  a  background  for  the  crystals.  A  one 
inch  strip  of  11  different  ranges  of  crystal,  with  approximately  1/2 
inch  of  black  paint  between  them,  was  applied  over  the  black  paint. 
Two  coats  of  each  crystal  were  applied,  using  a  small  paint  brush.  The 
first  coat  was  allowed  to  dry  completely  before  the  second  coat  was 
applied.  After  the  crystals  were  dry,  two  coats  of  Rez  polyurethane 
(gloss  clear  plastic  coating,  interior-exterior  77-5)  coating  were 
applied  directly  over  the  crystals.  The  polyurethane  coating  was 
applied  to  protect  the  crystals  from  wind  blown  sand  and  from  the 
ultraviolet  rays  of  the  sun.  Complete  liquid  crystal  calibration 
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results  are  shown  in  Table  1.  Details  of  the  liquid  crystal  research 
will  not  be  dwelt  upon  further  in  this  report.  The  interested  reader 
is  referred  to  Reference  2  for  additional  information. 

The  rocket  motor  storage  container  system  was  then  moved  to  the 
China  Lake  dump  storage  site.  The  system  was  aligned  in  a  north-south 
direction,  well  away  from  the  influence  of  other  ordnance  (Figure  6). 
The  thermocouple  leads  were  connected  through  a  junction  box  and  under- 
ground cable  to  a  Honeywell  Electronik  25  recorder  which  had  been 
calibrated  to  read  the  thermocouple  output  directly  in  degrees  Fahren- 
heit to  an  accuracy  of  ±  1°F.  The  recorder  was  located  in  an  air- 
conditioned  shed  about  60  feet  from  the  system. 

Initial  data  indicated  that  thermocouple  #7  was  not  responding 
properly,  and,  therefore,  its  output  was  neglected.  Initial  color 
photographs  were  taken  of  the  liquid  crystals,  and  it  was  immediately 
apparent  that  good  thermal  mappings  could  be  obtained  if  the  crystals 
were  stable  under  the  adverse  desert  environment.  The  brilliance  of 
the  colors  exhibited  by  the  crystals  under  the  bright  desert  sun  was 
much  better  than  had  been  expected.  The  system  was  allowed  two  weeks 
to  reach  a  periodic  steady  state  before  additional  photographic  data 
were  obtained. 

Extensive  photographic  data  were  collected  on  27  and  28  July  1972 
after  two  weeks  of  exposure  to  the  desert  environment.  Both  super  8  mm 
and  16  mm  color  movies  and  35  mm  color  slides  were  taken  of  the  liquid 
crystals.  No  colored  filters  were  used  on  any  of  the  cameras,  although 
standard  haze  filters  were  used  to  take  the  super  8  mm  movies  and  most 
of  the  35  mm  slides. 

At  this  time  a  second  storage  container,  this  one  without  a  rocket 
motor  inside,  was  instrumented  with  intrinsic  thermocouples  in  the  same 
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TABLE  I 
Calibration  of  Liquid  Crystals 

Calibration  Bath 
NCR                        Manufacturer's  2  Coats 

Desig.   Color  Change Responses Liquid  Crystals 


R-27 


R-33 


R-37 


R-41 


R-45 


R-49 


R-53 


R-56 


R-59 


S-62 


S-64 


Red 

Green 

Blue 

Red 

Green 

Blue 

Red 

Green 

Blue 

Red 

Green 

Blue 

Red 

Green 

Blue 

Red 

Green 

Blue 

Red 

Green 

Blue 

Red 

Green 

Blue 

Red 

Green 

Blue 

Red 

Green 

Blue 

Red 

Green 

Blue 

27.0  25.6+. 5 

28.6  28.0+.5 

30.0  28.7+. 5 

33.0  32.7+. 5 

34.6  33.3+. 5 

36.0  34.2+. 5 

37.0  36.2+. 5 

38.6  37.1+.5 

40.0  38.0+.5 

41.0  40.3+.5 

42.6  41.0+.5 

44.0  42.0+.5 

45.0  42.8+. 5 

46.6  43.6+. 5 

48.0  44.3+. 5 

49.0  46.7+. 5 

50.6  47.1+.5 

52.0  48.4+. 5 

53.0  50.5+.5 

54.6  52.1+.5 

56.0  53.3+. 5 

56.0  53.8+. 5 

57.6  56.0+.5 

59.0  56.5+. 5 

59.0  56.9+. 5 

60.6  57.5+. 5 

62.0  58.9+. 5 

62.0  60.1+.5 

62.6  60.4+.5 

63.0  60.9+.5 

64.0  60.9+.5 

64.6  61.4+.5 

65.0  62.7+. 5 
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manner  as  the  previous  container.  As  only  three  data  channels  remained 
open  on  the  recorder,  only  three  thermocouples  were  applied  to  this  new 
container.  The  three  thermocouples  were  applied  at  the  0300,  0900,  and 
1200  positions  at  the  axial  midpoint  of  the  container.  This  container 
was  set  end  to  end  with  the  system  that  was  already  in  place  at  the  site. 
The  purpose  of  this  study  was  to  determine  if  the  inclusion  of  the  rocket 
motor  in  the  container  had  a  significant  effect  on  the  surface  tempera- 
ture of  the  container.  Thermocouple  #7  was  connected  at  the  0900  posi- 
tion, #23  at  the  1200  position,  and  #24  at  the  0300  position.  It  was 
immediately  apparent  that  thermocouple  #7  was  continuing  to  give 
unreliable  readings,  and  therefore,  the  data  taken  on  channel  #7  were 
again  neglected. 

IV.  EXPERIMENTAL  RESULTS 

The  data  presented  in  Tables  2  and  3  were  obtained  from  the  thermo- 
couples on  the  rocket  motor  storage  container  system  located  at  China 
Lake,  California.  The  thermocouple  output  was  read  out  on  a  Honeywell 
Electronik  25,  24  channel  recorder  which  had  been  calibrated  at  50,  100, 
and  150°F.  The  data  were  taken  on  two  consecutive,  typical  summer  days 
(August  1  and  2,  1972)  at  China  Lake.  (Figure  7  shows  the  location  of 
the  thermocouples  used  to  collect  this  temperature  data.)  Each  thermo- 
couple was  read  once  every  24  minutes.  The  data  shown  in  Table  2 
present  the  storage  container  temperature  at  four  locations  plus  three 
different  ways  of  averaging  this  data.  It  also  presents  the  ambient 
temperature  and  the  approximate  time  of  day.  The  data  shown  in  Table  3 
present  the  surface  temperature  of  the  rocket  motor  and  three  ways  to 
average  this  data.  Also  shown  are  the  temperature  at  the  center  of  the 
rocket  motor  and  the  approximate  time  of  day. 
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Figure   7:    Thermocouple  Locations  for  Experimental  Data 
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Table  4  compares  the  temperatures  measured  at  the  12  o'clock  loca- 
tion on  the  empty  container  with  those  measured  on  the  container  which 
contained  the  ASROC  motor.  For  convenience,  in  the  same  table,  the 
bulk  container  temperature  has  been  compared  with  the  average  of  the 
12  o'clock  temperature  and  the  ambient  temperature.  The  reason  for 
presenting  this  information  is  discussed  under  Conclusions  and  Recom- 
mendations. 

V.  THEORETICAL  ANALYSIS 

A.  Introductory  Remarks 

In  the  present  analytical  investigation  the  container  surface 
temperature  variation  was  assumed  known  a  priori.  The  results  of  our 
experimental  investigation  indicate  that  there  is  a  significant  tempera- 
ture variation,  as  much  as  45° F,  around  the  circumference  of  the 
container.  The  heat  transferred  to  the  rocket  motor  stored  inside  the 
container  arrives  by  the  combined  effects  of  thermal  radiation  and 
convection  across  the  air  gap  separating  the  motor  and  the  container. 
The  air  gap  represents  a  significant  thermal  resistance  to  the  heat 
transfer  process  and,  further,  tends  to  smooth  the  spatial  temperature 
variation  around  the  circumference  of  the  motor.  This  may  be  seen  by 
examining  the  experimental  data,  Table  4,  and  noting  that  the  motor 
circumferential  temperature  variation  is  a  maximum  of  18°F,  and  this 
only  occurs  for  a  short  time  span  each  day.  Accordingly,  it  seemed 
reasonable  to  assume  that  within  the  rocket  motor,  heat  transfer  would 
be  essentially  one-dimensional,  with  the  primary  direction  being  radial. 
In  order  to  develop  a  one-dimensional  analytical  solution,  the  container 
circumferential  temperature  variation  was  averaged  to  develop  a 
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TABLE  4 

Comparison  of  Measured  Temperatures  on  two  ASROC  containers: 
One  container  contained  a  sand-filled  motor,  and  the  second 

container  was  empty 


Time 

Ambient 

Top, 
Filled  Container 

Top, 
Empty  Container 

Average 
Air&Contain. 

Actual 
Bulk 

Aug. 1,1972 

(°F) 

(°F) 

(°F) 

(°F) 

(°F) 

0536 

76 

69 

67 

73 

75 

0600 

77 

68 

66 

73 

75 

0624 

80 

73 

70 

77 

78 

0648 

80 

81 

79 

81 

83 

0712 

83 

90 

88 

87 

88 

0736 

85 

98 

98 

92 

92 

0800 

87 

105 

105 

96 

96 

89 

110 

no 

99 

99 

91 

116 

116 

104 

102 

92 

121 

121 

107 

106 

94 

129 

128 

111 

110 

1000 

96 

133 

133 

114 

113 

97 

139 

139 

118 

116 

100 

143 

143 

121 

119 

101 

147 

146 

124 

122 

103 

150 

147 

126 

124 

1200 

104 

153 

150 

128 

126 

106 

156 

152 

131 

130 

106 

157 

152 

131 

131 

106 

154 

150 

130 

132 

107 

157 

150 

132 

134 

1400 

108 

154 

150 

131 

136 

109 

153 

149 

131 

136 

no 

148 

145 

129 

134 

110 

143 

138 

127 

133 

109 

142 

137 

126 

133 
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TABLE  4  (Continued) 

Top,         Top,        Average    Actual 
Time      Ambient  Filled  Container  Empty  Container  Air&Contain.  Bulk 


Aug. 1,1972) 
cont'd.) 

1600 

107 

137 

133 

122 

132 

108 

134 

130 

121 

129 

108 

130 

126 

119 

128 

106 

126 

121 

116 

126 

• 

104 

123 

120 

114 

125 

1800 

103 

118 

114 

111 

121 

101 

113 

109 

107 

118 

99 

107 

104 

103 

113 

96 

101 

99 

98 

105 

95 

91 

88 

93 

96 

2000 

93 

88 

84 

91 

94 

91 

86 

82 

89 

92 

90 

85 

82 

88 

91 

89 

84 

80 

87 

90 

87 

83 

79 

85 

88 

2200 

87 

81 

78 

84 

86 

85 

80 

77 

82 

85 

86 

79 

76 

82 

84 

86 

79 

76 

82 

84 

84 

79 

76 

81 

84 

2  Aug 
0000 

82 

77 

74 

79 

82 

81 

74 

70 

77 

80 

81 

72 

67 

76 

79 

81 

72 

67 

76 

79 

79 

72 

68 

75 

78 

0200 

79 

72 

68 

75 

78 

79 

72 

68 

75 

77 

78 

72 

68 

75 

77 

77 

71 

67 

74 

76 

75 

70 

66 

72 

75 
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TABLE  4  (Continued) 

Top,         Top,        Average   Actual 
Time      Ambient  Filled  Container  Empty  Container  Air&Contain.  Bulk 


Aug. 2, 1972) 
cont'd.) 

0400 

73 

67 

63 

70 

73 

72 

66 

61 

69 

72 

73 

65 

61 

69 

71 

69 

64 

60 

67 

70 

69 

63 

58 

66 

69 

0600 

68 

63 

58 

65 

68 

71 

65 

59 

68 

71 

77 

74 

70 

76 

76 

79 

84 

83 

82 

82 

80 

93 

92 

87 

88 

0800 

83 

101 

99 

91 

92 

85 

108 

108 

96 

96 

86 

114 

114 

100 

100 

88 

121 

121 

104 

104 

89 

126 

126 

107 

106 

1000 

91 

130 

129 

109 

108 

94 

137 

135 

115 

113 

96 

143 

141 

119 

116 

100 

146 

144 

123 

119 

100 

149 

147 

124 

122 

1200 

101 

151 

148 

126 

124 

103 

150 

147 

127 

126 

104 

156 

150 

130 

129 

104 

154 

149 

129 

130 

105 

155 

147 

130 

132 

1400 

105 

156 

149 

131 

134 

107 

149 

142 

128 

133 

107 

149 

142 

128 

134 

106 

151 

143 

129 

138 

107 

144 

138 

126 

137 

1600 

108 

139 

134 

124 

135 

107 

136 

132 

122 

133 
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TABLE  4  (Continued) 

Top,  Top,       Average   Actual 

Time      Ambient  Filled  Container  Empty  Container  Air&Contain.  Bulk 


Aug. 2, 1972) 
cont'd.) 

107 

136 

130 

122 

133 

105 

131 

126 

116 

132 

103 

124 

118 

114 

129 

1800 

103 

117 

113 

110 

123 

100 

111 

109 

106 

116 

- 

98 

105 

102 

102 

112 

95 

99 

96 

97 

103 

93 

91 

88 

92 

96 

2000 

91 

87 

84 

89 

93 

90 

85 

82 

87 

91 

88 

83 

80 

86 

89 

87 

81 

77 

84 

87 

85 

80 

75 

82 

86 

2200 

84 

78 

72 

81 

84 

82 

75 

71 

79 

82 

85 

74 

67 

79 

81 

83 

76 

70 

79 

81 

81 

74 

69 

77 

80 

3  Aug. 
0000 

78 

71 

66 

75 

78 

81 

69 

62 

75 

77 

79 

70 

66 

75 

77 

79 

72 

70 

75 

77 

78 

71 

67 

74 

76 

0200 

77 

70 

66 

73 

75 

73 

69 

66 

71 

74 

72 

65 

60 

69 

72 

70 

64 

59 

67 

71 

67 

63 

57 

65 

69 

0400 

67 

61 

56 

64 

68 

66 

60 

55 

63 

67 

65 

60 

54 

63 

66 
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TABLE  4  (Continued) 

Top,  Top,        Average   Actual 

Time      Ambient  Filled  Container  Empty  Container  Air&Contain.  Bulk 


Aug. 3,1972) 
cont'd.) 

65 

58 

52 

61 

65 

65 

58 

51 

61 

64 

0600 

65 

58 

52 

61 

64 
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relationship  between  the  container  bulk  temperature  and  time.  By 
definition: 


Tbulk    r 


cpTdV 


(1) 


cpdV 


where 


T   =   container  temperature 
c   =   container  specific  heat 
dV   =    container  elemental  volume 


Noting  that  the  specific  heat  of  steel,  c  ,  is  essentially  constant 
over  the  temperature  range  of  interest  and  that  the  elemental  volume, 
dV,  is  directly  proportional  to  the  elemental  angular  increment,  dy, 
Equation  (1)  reduces  to: 


Tbulk  ■  XTJ         Td*  <2> 

Since  container  surface  temperature  data  were  available  only  at 
discrete  locations  around  the  circumference,  equation  (2)  was  approxi- 
mated as: 

Tbulk  *  ?7-LTi(i-)  (3) 

i=l 


i=l 
where  the  T.'s  are  the  measured  temperatures  at  the  12,  3,  6,  and  9 

o'clock  positions  on  the  container. 
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Equation  (3)  was  used  to  convert  local  experimental  data  to  bulk 
temperature  information.  Figure  8  shows  the  bulk  temperature  variation 
as  a  function  of  time  as  experimentally  recorded  on  2  August  1972. 
Shown  as  a  continuous  curve  on  the  same  graph  is  a  pure  sine  wave  with 
the  same  period  and  amplitude.  While  the  agreement  is  not  perfect, 
note  that  over  a  significant  portion  of  the  transient,  the  agreement  is 
quite  close.  Due  to  this  agreement  it  was  decided  to  model  the  container 
bulk  temperature- time  variation  as  a  pure  sine  wave  of  the  form: 

Tc  *  (Tm  "  V  sin  ut  +  Ta  <4> 

where 

T   =  container  bulk  temperature 

T   =  maximum  bulk  temperature 

T   =  average  bulk  temperature 
a 

t   =  time 

a)   =  circular  frequency  (2u/period) 

B.  One-dimensional  Analytical  Model 

The  method  of  complex  temperatures,  as  presented  by  Arpaci  [1], 
was  used  to  find  the  radial  temperature  distribution  in  a  cylindrical 
rocket  motor  stored  in  a  cylindrical  container.  The  container  tempera- 
ture was  assumed  to  vary  in  the  form  given  by  Equation  (4).  The  complete 
analytical  derivation  may  be  found  in  reference  2.  The  major  assumptions 
employed  in  the  model  are:  heat  flow  is  one  dimensional  in  the  radial 
direction  only,  no  heat  sources  or  sinks  exist  within  the  motor,  and 
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the  sinusoidal ly  varying  surface  temperature  is  spacially  uniform  over 
the  entire  container  surface. 

Heat  transfer  across  the  air  gap  was  treated  by  employing  the 
concept  of  an  overall  gap  conductance,  h.  This  overall  conductance 
combines  the  effects  of  both  radiation  and  convection.  Assuming  that 
the  radiation  and  convection  processes  are  uncoupled,  one  may  write: 


h  ■  hr  +  hc  (5) 


where 


h  =  linearized  gap  radiation  coefficient 
h  ■  gap  convective  coefficient 
The  net  heat  transfer  rate  between  the  motor  and  container  is: 

Ql-2   s   *A1  (T1  "  V  (6) 


where 


Q,  2   -   heat  transfer  rate  between  motor  and 
container 


A,    =   motor  surface  area 

T,    n      motor  bulk  surface  temperature 

T«    s   container  bulk  surface  temperature 

The  gap  radiation  coefficient,  h  ,  was  found  by  linearizing  the  grey 
body  radiation  relationship: 


33 


V2    =        Fl"2  °  Al   (T1     '  T2  > 


where 


hr 


Fl-2°   (T1  +T2)   (T12  +  T22)  Al    (Tr"  V         (7) 


hrAl   (T1  "  T2) 


Fl-2  °  (T1  +  V   (T12  +  T22)  (8) 


and 


F,_2  =  Grey  body  exchange  factor 

o  =  Stefan-Boltzmann  constant 

T,  ■  Motor  skin  temperature 

T2  =  Container  temperature 

The  gap  convective  coefficient,  h  ,  may  be  found  by  using  the 
correlations  developed  by  Liu,  Mueller,  and  Landis  [3]  for  convection 
across  cylindrical  gaps,  namely: 


where 


(0.278) 

hc  -       [°-135kq   ]        [Pr2fir] 


k  =        Thermal  conductivity  of  the  gas  in  the  gap 

Pr  =        Prandtl  number  =  v/a 

v  =        Gas  kinematic  viscosity 
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(9) 


Gas  thermal  diffusivity 


Gr   =   Grashof  number  =  g3AIAr 


g   =   Local  gravity  acceleration 


6   =   Thermal  coefficient  of  expansion 


AT   =    (T,  -  T2) 
Ar   =    (r2  -  r}) 


r,   ■   motor  radius 

r«   s   container  radius 

The  steady,  periodic  solution  describing  the  radial  temperature 
distribution  in  the  rocket  motor,  T(r,t),  is  developed  in  detail  in 
Reference  2  using  the  method  of  complex  temperatures  [1],  The  solution, 
in  normalized  form,  is: 


(T(r,t)  -  Ta) 
<Tm  -  V 


[be/ (As)  +  bei^(Ac) 
rx^  +  X.^jl/2 


1/2 


sin  (wt  +  6)  (10) 


or 


e(r,  t) 


e  sin  (ut  +  6) 


(11) 


where 


T(r,t) 

Ta 
Tm 


Motor  temperature  at  radial  location  r  and  time  t 


Average  container  temperature 


Maximum  container  temperature 
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far.t)  -  T 
e(r,t)  =        Normalized  temperature    -p-= = =- 


L  m 

t  -        Time 


Radial  location 


r        -        Motor  radius 


Normalized  radial  location  (r/r  ) 


u)        ■    Circular  frequency  of  container  temperature 
variation 


V? 


A  s  Normalized  frequency 

Thr  ] 

B  -  Biot  number        -r~-\ 

k  =  Motor  thermal  conductivity 

p  =  Motor  density 

c  =  Motor  specific  heat 

a  =  Motor  thermal  diffusivity  (k/pc) 

ber  =  Real  Bessel  function 

bei  s  Imaginary  Bessel  function 

X„  *  ber  (A)     +    — *—    ber, (A)     +    -A_    bei, (A) 
r  °  B\/2  B>/2 

Xi  =  beiQ(A)     +        A        bei^A)     -        A        ber^A) 

B\/2  B\/2 
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.1 
-    tan 


r      o 

(Ac) 

"Xi 

ber0 

(AC)" 

KbQrn 
r      o 

(Ad 

+  Xi 

bei 

0 

(AC) 

r  2        2  i 

e  ■  Temperature  amplitude     ber  (Ac)  +  bei  (Ac) 
attenuation  L  r  2 ^T? 

[Xr  +  Xi] 


/2 


A  parameter  study  was  carried  out  on  equation  (10)  with  the 
assistance  of  an  IBM  360  computer.  Values  of  the  temperature 
amplitude  attenuation  parameter,  e  ,  and  phase  lag,  6,  were 
determined  at  three  normalized  radial  locations,  c  =  0.0,  0.5, 
and  1.0,  for  values  of  A  and  B  of  practical  interest.  Tables  5 
through  9  present  the  results  of  this  study.  The  phase  lag,  6, 
is  given  in  radians,  where  2ir  radians  equals  one  complete  cycle. 
For  environmental  considerations,  one  complete  cycle  represents 
twenty- four  hours. 

For  ease  of  use,  the  results  of  the  parameter  study  are  also 
presented  in  nomograph  form  in  Figures  9  through  12.  Temperature 
amplitude  attenuation,  e  ,  and  phase  lag,  6,  are  plotted  parametri cal ly 
against  A  and  B  at  normalized  radial  locations  of  c  =  0.0  and  c  x  1.0 
(center  and  surface  of  motor,  respectively). 
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Figure  12.  Variation  in  Phase  Lag,  6  at  the  Motor  Center,  with 
Biot  number,  B,  for  various  values  of  the  Normalized 
Frequency  Parameter,  A. 
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C.  Sample  Calculation 

As  an  example  of  the  use  of  the  temperature  attenuation  and  phase 
lag  nomographs,  consideration  will  be  given  the  once-fired  ASROC  motor 
used  in  the  experimental  investigation.  The  motor  was  filled  with  dry, 
wind  blown  sand.  Reference  4  lists  the  properties  of  sand  as: 

p  *  94.8  lbm/ft3 

k  =  0.188  BTU/hr-ft-°F 

c  =  0.195  BTU/lbm-°F 


a  =  0.01  ft  /hr 


The  thermal  properties  of  the  air  in  the  gap  between  the  motor 
and  container  were  taken  from  reference  5.  An  average  air  temperature 
of  100°F  was  used  as  a  representative  reference  temperature. 

Pr  =  0.72 

6  3.  -1 


£§■«  1.76  x  10  (°F  -  ft  ) 

k    =    0.0154  BTU/hr-ft-°F 

The  gap  convective  coefficient  defined  by  equation  (9)  was  calcu- 
lated using  an  estimated  temperature  difference  between  motor  and 
container  of  10°F.     The  results  of  the  calculation  yield: 

h  n  ?n  BTU 

The  linearized  gap  radiation  coefficient,  equation  (8),  was 
calculated  using  the  following  values: 
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r,       =     0.5  ft. 


r2      =    0.75  ft. 


! 


o         =     0.1714  x  10"8  hr.ft".efti, 

T,       =     100°F  ■  560° R    (average  motor  temperature) 

T2      =     100°F  =  560° R    (average  container  temperature) 

e1  =  e2  -  0.9  (typical  emissivities  for  grey  painted  surfaces) 

F,  0  =    J =     0.84  (ref.   5,  p.  260) 

1-2         _UlW_L-i\ 
el       r2    \e2        / 

Substituting  these  values  into  equation  (8),  one  finds 

h    -  ini;       BTU 

Using  the  calculated  values  of  h  and  h  ,  the  overall  gap 
conductance  is  found  to  be: 


h    =  hr+hc 


i  ?k  BTU 


The  frequency  parameter,  A,  and  Biot  number,  B,  are  now  easily 
calculated  as  A  =  2.60  and  B  =  3.35.     In  calculating  A,  the  circular 
frequency  was  assumed  to  equal  2*/24  =  0.262    rjj^ns     . 

Now,  referring  to  Figures  9  through  12,  the  following  values  are 
found : 
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%    =  0.0,  ep  =  0.42  6  =  -1.8  radians  (12) 

C  s  1.0,  er  *  0.67  6  =  -0.39  radians  *       (13) 

Substituting  the  above  values  into  equation  (11)  we  find: 
At  c  =  0  (motor  center) 


(T(o,t)  -  T  ) 


^—  =  0.42  sin  (at  -  1.8)  (14) 


<Tm  -  V 
At  c  ■  1.0  (motor  skin) 

(T(ro»t}  *  V  (U) 

iT°  :  T  v  a   *  0.67  sin  (u.t  -  0.39)  Ub; 

To  convert  e  and  6  to  dimensional  temperature  and  time,  we  must 
specify  the  following  information: 

1.  The  duration  of  one  complete  cycle. 

2.  The  maximum  container  bulk  temperature. 

3.  The  minimum  container  bulk  temperature. 

4.  The  average  container  bulk  temperature. 

As  an  example,  using  the  experimental  results  obtained  on  August  2, 
1972  (Table  2),  we  find  that  one  complete  cycle  is  equivalent  to 
twenty-four  hours  and  the  maximum  and  minimum  container  bulk  tempera- 
tures are  138°F  and  64°F,  respectively.  These  temperatures,  when 
averaged,  yield  an  average  container  temperature  of  101 °F.  Substituting 
these  results  into  equations  (14)  and  (15)  and  rearranging  in  dimensional 
form,  one  finds: 

At  r  =  0.0  ft.  (center  of  motor), 

T(o,t)  =  15.5  sin  {(0.26)  (t  -  6.9)}  +  101  (°F)  (16) 
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At  r  =  0.5  ft  (surface  of  motor), 

T(rQ,t)  =  25  sin  {(0.26)  (t-1.5)}  +  101  (°F)  (17) 

Figure  13  shows  plots  of  equations  (16)  and  (17)  together  with  the 
assumed  container  bulk  temperature  variation.  The  important  points  to 
note  are: 

1.  Maximum  and  minimum  container  bulk  temperatures 

of  138°F  and  64°F,  respectively,  lead  to  predicted 
maximum  and  minimum  temperatures  at  the  motor  center 
of  116. 5° F  and  84. 5° F,  respectively,  and  at  the  motor 
surface  of  126°F  and  76° F,  respectively. 

2.  The  thermal  time  lag  at  the  rocket  motor  center  is 
predicted  to  be  6.9  hours  while  that  at  the  motor  surface 
is  1.5  hours. 

VI.  COMPARISON  OF  THEORY  WITH  EXPERIMENT 

Figures  14  and  15  show  comparisons  of  theoretically  predicted 
temperatures  with  experimental  results  obtained  on  August  2,  1972. 
The  experimentally  determined  container  bulk  temperature  variation 
was  used  as  the  forcing  function  in  the  theory.  Figure  14  represents 
the  bulk  temperature  variation  at  the  surface  of  the  motor  while 
Figure  15  represents  the  temperature  variation  at  the  center  of  the 
motor. 

Table  10  compares  theoretical  amplitude  and  phase  shift  predictions 
at  the  motor  center  and  skin  with  experimental  results.  As  can  be  seen, 
the  comparison  is  excellent  with  the  exception  of  the  phase  shift  com- 
parisons at  the  motor  center.  An  uncertainty  analysis  has  been  performed 
on  both  the  experiment  and  theory  (the  uncertainty  in  the  theory  is  a 
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Figure  13.  Analytical  Predictions  of  Temperature  Variation  with 
Time  at  surface  and  center  of  ASR0C  motor. 
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Figure  14.  Comparison  of  Measured  Values  of  Motor  Skin  Temperatures 
with  Values  Predicted  Using  Equation  17. 
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Figure  15.  Comparison  of  Measured  Values  of  Motor  Center  Temperatures 
with  Values  Predicted  Using  Equation  16. 
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result  of  uncertainties  in  the  values  of  the  thermophysical  properties 
of  sand  used  in  the  calculations),  and  it  is  estimated  that  the  total 
uncertainty  may  be  represented  with  uncertainty  bands  of -±3°F  and  ±0.5 
hours  on  each  data  point  (see  reference  2  for  details).  Theory  and 
experiment  agree  within  these  uncertainty  limits  at  values  of  maximum 
and  minimum  temperature. 

Table  10 

Comparison  of  Experimentally  Determined  and 
Theoretically  Predicted  Values  of  Temperature  and 
Phase  Shift  on  the  Surface  and  in  the  Center  of  a 
Sand-filled  ASROC  motor 


Radial 
Location 

Maximum  Bulk 
Temperature  (°F) 

Minimum  Bulk 
Temperature  °F 

Phase  Shift 
(hours) 

Theory 

Experiment 

Theory 

Experiment 

Theory 

Experiment 

Motor 
Center 

116.5 

117 

84.5 

84 

6.9 

5.3 

Motor 
Skin 

126 

123 

76 

75 

1.5 

1.6 

The  relatively  poor  agreement  between  theory  and  experiment  during 
certain  portions  of  each  cycle  may  be  attributed  to  several  shortcomings 
of  the  analytical  model: 

1.  The  model  assumes  one  dimensional  heat  flow.  In  reality, 
the  flow  is  two  dimensional  and  possibly  mildly  three 
dimensional. 

2.  The  model  assumes  that  the  gap  conductance  is  constant.  In 
reality,  the  gap  conductance  varies  with  time  due  to  the 
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fact  that  the  temperature  difference  between  the  motor  skin 
and  container  varies  with  time. 

3.  The  model  assumes  that  the  bulk  container  temperature  varies 
si nusoi dally  with  time.  In  reality,  the  variation  is  not 
truly  sinusoidal  as  can  be  seen  by  examining  Figure  8. 

4.  In  order  to  develop  a  one-dimensional  solution,  a  bulk  container 
temperature,  Equation  (3),  was  defined  and  used  as  the  forcing 
function  in  the  theoretical  model.  This  bulk  temperature  only 
partially  compensates  for  the  truly  two  dimensional  nature  of 
the  problem. 

All  of  the  aforementioned  shortcomings  of  the  one-dimensional 
analytical  model  can  be  overcome  by  employing  a  numerical  procedure  to 
solve  for  the  temperature  distribution  in  the  ordnance.  This  has  been 
done  and  is  reported  on  in  detail  in  reference  2.  While  the  numerical 
solution  yields  results  that  very  closely  approximate  those  obtained 
experimentally,  the  increase  in  time  and  effort  is  considerable  when 
compared  to  the  simplicity  of  using  the  analytical  model,  especially 
in  nomograph  form. 
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VII.  CONCLUSIONS  AND  RECOMMENDATIONS 

Based  on  the  results  of  our  investigation  we  offer  the  following 
conclusions  and  recommendations: 

1.  Although  a  sine  wave  is  not  a  perfect  fit  to  the  actual 
container  temperature  variation  at  all  points  in  time,  it  is  useful 
for  predicting  bulk  temperatures  and  thermal  time  lags  throughout 
the  rocket  motor,  especially  if  only  high  and  low  bulk  temperatures 
are  of  interest. 

2.  The  simplicity  of  the  analytical  model,  as  presented  in  nomo- 
graph form,  allows  parameter  studies  to  be  quickly  and  easily  carried 
out.  The  results  of  our  study  should  prove  to  be  a  useful  tool  to  aide 
in  the  future  design  of  thermally  optimized  containers. 

3.  Our  theoretical  model  has  been  verified  by  comparison  with 
experimental  results  obtained  on  a  container  stored  ASROC  motor. 
Additional  comparisons  should  be  made  with  data  obtained  on  other  types 
of  container  stored  ordnance.  Vast  amounts  of  experimental  data  have 
been  obtained  on  a  variety  of  ordnance  and  are  available  for  compari- 
son [6].  However,  before  our  model  can  be  employed,  a  method  for 
obtaining  meaningful  container  bulk  temperature  variations  must  be 
developed.  Virtually  all  existing  data,  with  the  exception  of  that 
presented  herein,  were  taken  at  the  12  o'clock  locations  on  the  container 
and  motor.  While  these  temperatures  represent  extremes,  or  near  extremes, 
they  are  not  truly  representative  of  bulk  temperatures.  A  technique  is 
needed  for  estimating  maximum  and  minimum  bulk  temperatures  based  on  a 
knowledge  of  only  the  temperatures  taken  at  the  12  o'clock  locations 
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and  the  measured  ambient  temperature,  the  measured  ambient  temperature 
being  somewhat  representative  of  the  temperature  at  the  6  o'clock 
location  on  the  container.  As  an  example,  Table  4  compares  the  actual 
bulk  temperature  of  the  ASROC  container  used  in  the  present  studies 
with  an  estimated  bulk  temperature  calculated  by  simply  averaging  the 
temperature  measured  at  the  12  o'clock  location  with  the  ambient  air 
temperature.  We  do  not  necessarily  recommend  this  scheme  in  prefer- 
ence to  other  schemes  but  simply  offer  it  as  an  example.  Perhaps  a 
weighted  averaging  technique  would  yield  closer  agreement  as  the 
straight  averaging  technique  does  not  yield  close  comparisons  at  all 
points  in  time,  especially  near  peak  temperatures. 

4.  Table  4  presents  a  comparison  of  the  temperatures  measured  at 
the  12  o'clock  locations  on  the  container  which  contained  an  ASROC 
motor  with  a  container  that  was  empty.  Close  inspection  of  the  data 
indicate  that  there  is  very  little  difference  in  the  temperatures 
measured  on  the  two.  This  suggests  that  an  empty  storage  container 

may  be  used  to  obtain  container  temperature  data  for  use  in  the  analyti- 
cal model. 

5.  The  analytical  model  presented  herein  uses,  as  a  forcing  function, 
experimentally  determined  values  of  the  container  temperature.  It  is 
recommended  that  a  study  be  conducted  to  develop  predictions  of  the 
container  surface  temperature  as  a  function  of  solar  irradiation, 
terrestial  irradiation,  wind  velocity,  and  ambient  air  temperature. 
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